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We show that an inhomogeneous cosmological lepton number may have produced turbulence in 
the primordial plasma when neutrinos entered (almost) free-streaming regime. This efltect may 
be responsible for the origin of cosmic magnetic fields and give rise to a detectable background of 
gravitational waves. An existence of inhomogeneous lepton asymmetry could be naturally generated 
by active-sterile neutrino oscillations or by some versions of Affleck and Dine baryogenesis scenario. 



The idea that the early universe went through one or 
more turbulent phases is recurrent in the scientific lit- 
erature since the beginning of the modern cosmology. 
Among other reasons, cosmic turbulence was often in- 
voked to explain the origin of the magnetic fields (MFs) 
observed in most spiral galaxies and galaxy clusters . 
One of the first attempts Q was based on the observa- 
tion that in the radiation era a weak MF should be gener- 
ated by turbulent eddies because the rotational velocities 
of relativistic electrons and non-relativistic ions would 
change differently during universe expansion. Subse- 
quently, turbulence may have induced a Magneto-Hydro- 
Dynamical (MHD) dynamo which amplified exponen- 
tially the field until equipartition between the plasma 
turbulent kinetic energy and the MF energy is eventually 
reached. This nice scenario, however, remains unaccom- 
plished in the absence of a plausible mechanism of tur- 
bulence generation. Since gravitational forces conserve 
angular moment, vorticity must have been produced by 
non-gravitational forces. Primordial phase transitions 
might have produced MFs if they were first order [Q. 
The idea is that the expanding bubble walls between two 
phases give rise to small electric currents which power 
a seed MF. Turbulence appears when bubbles collide at 
the end of the phase transition, or hydrodynamic insta- 
bilities develop on the bubble walls, producing MHD am- 
plification of the seed field. Several applications of this 
idea have been studied for the case of the quark-hadron 
(QHPT) and electroweak (EWPT) phase-transitions (see 
1^ and refs. therein). Unfortunately, it is still unclear if 
any of these transitions can really be first order. Fur- 
thermore, the major problem with this kind of scenario, 
is that it can hardly account for large scale MFs. The 
reason is that the comoving coherence length of the MF 
is at most given by the Hubble horizon at the phase 
transition which is much smaller than a typical galaxy 
size. Although the coherence length may grow due to 
MHD effects this generally happens to the expenses of 
the MF strength. Careful studies showed that the EWPT 



and QHPT cannot account for the galactic and cluster 
MFs To overcome the small scale problem several 
mechanisms of generation of seed magnetic fields at in- 
flationary stage have been proposed which are operative 
if conformal invariance of electromagnetic interaction is 
broken (for more references see the review |2]). 

In this Letter we propose a new mechanism for the 
generation of cosmic MFs. In contrast to many previ- 
ously considered mechanisms it does not demand any 
new strong physical assumptions and can be realized, in 
particular, in almost minimal standard model of particle 
physics. It might operate at neutrino decoupling epoch, 
T ~ 1 — 2 MeV, i.e. when the Hubble horizon was consid- 
erably larger than what it was at EWPT and at QHPT. 
Our basic assumption is that the net lepton number den- 
sity (iVa(x) = ni,^(x) — np^(x), a = e,/i, r) of one, or 
more, neutrino species was not uniform before neutrino 
decoupling and changed in space over some characteristic 
scale A which could be smaller than the Hubble horizon 
at the decoupling time. 

As a result, when the neutrino mean free path (-v{T) 
grew and became comparable to A, neutrino currents 
should be developed along the density gradients. We will 
show that elastic scattering of the diffusing neutrinos on 
electrons and positrons would be able to accelerate the 
electron-photon fluid producing vorticity in the plasma. 
Turbulence may have developed by this process in the 
short interval of time during which the random forces due 
to neutrino elastic scattering overcome the shear viscos- 
ity force. Depending on the amplitude and wavelength 
of the fluctuations of A^a(x), this period could be suf- 
ficient for the MHD engine to generate magnetic field 
in equipartition with the turbulent kinetic energy. The 
seed field required to initiate the process arises naturally 
as a consequence of the difference between the t'aC" and 
VaG^ cross sections and of the neutrino-antineutrino local 
asymmetry. Furthermore, we will show that turbulence 
should give rise to a background of gravitational waves. 

The first step of our computation is to determine the 
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neutrino momentum flux produced by the inhomoge- 
neous -/Vq(x). Starting form the Boltzmann equation we 
obtain the equation describing evolution of the average 
flux of i-th component of neutrino momentum 
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where 
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In the above E and k are respectively the neutrino en- 
ergy and spatial momentum, H is the universe expan- 
sion rate, is the effective weak interaction time, and 
fv{E, k, X, t) is the neutrino distribution function. As we 
see in what follows, the source of the magnetic field is 
proportional to curl of electric current, V x J, which in 
turn is proportional to the local vorticity of the source 
term diKij in Eq.(0). The latter is nonvanishing for 
anisotropic random initial distribution of neutrino lep- 
tonic charge and is numerically close to Kij divided by 
A. Hence, Eq.(||) becomes 
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-4i//CA(x,t)-r-^/CA(x,i) , 



where /C(x, t) = K/pj^ is the specific momentum flux on 
the length-scale A, and F^" is a unit vector parallel to 
the initial value of the vector djKij, which has a non- 
zero vorticity, i.e. V x 7^ 0. The exponential in the 
first term of the r.h.s. of Eq.(^) has been inserted to ac- 
count for the damping of the fluctuations due to neutrino 
diffusion. The last term in the r.h.s. of Eq.(||) represents 
the force per unit mass between the drifting neutrinos 
and the electron-positron fluid due to their scattering. 
It is understood that A changes with time due to uni- 
verse expansion. Since at T ~ 1 MeV, the Compton 
scattering rate was much larger than the universe expan- 
sion rate H, electrons, positrons and photons formed a 
(almost) single relativistic fluid at neutrino decoupling 
time (the role of baryons was negligible at this epoch). 
We assume the electron-photon fluid was approximately 
homogeneous before neutrino decoupling. This is possi- 
ble if the local excess of neutrinos of a given species is 
balanced by a deficiency of neutrinos of a different species 
(also sterile). The Euler equation of the fluid is 
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where p and p are the energy density and pressure of the 
electron-photon fluid, 7 

collision time, and rj = 4p^£^/15(p -I- p) is the shear vis- 
cosity due to neutrino diffusion. This approximation for 



— v'^ is the Lorentz factor, 
+ is the neutrino-electron(positron) mean 



the viscosity is valid in the limit when the neutrino mean 
free path is smaller than the characteristic scale of the 
problem. In our case this is not always true and the vis- 
cosity can be smaller. In Fig.l the macroscopic velocity v 
of the electron-photon fluid (continuous line) is presented 
as a function of the ratio x = t/td {td is the neutrino de- 
coupling time). It is obtained by numerical solution of 
the coupled set of Eqs. (||) and (||) for a particular choice 
of Srii, and characteristic size of the fluctuations (see the 
figure caption). It is visible from this figure that the ve- 
locity grows rapidly when the neutrino mean free path 
becomes comparable to the fluctuations size A and it is 
suppressed soon later when the viscosity on that length 
scale becomes dominant. 
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FIG. 1. The local fluid velocity v (continuous line), the 
specific neutrino momentum flux /C (dashed line) , and the ra- 
tio h = \e\B/T^ (dotted-dashed line) as functions of the time 
parameter x. We assumed here bn^jni, — 1 over a comov- 
ing length scale = {\/ H~^)t=tj, = 10"'^. The quantities 
plotted here represent avarages on the same scale. 

In this Letter we do not investigate the velocity and 
MF power spectrum. We only observe here that due to 
the rapid grow of the viscosity with time, the turbulence 
spectrum may be quite different from the Kolmogoroff's 
one. 

We now turn to the problem of magnetic field genera- 
tion. Starting from Maxwell equations in the primordial 
plasma and neglecting terms proportional to the inverse 
of the electric conductivity we obtain after simple algebra 



dtB + 2i/B = V X (v X B) + k^^V x Jext 



(4) 



where k ^ T/a |^,|| is the electric conductivity of the 
relativistic cosmological plasma {a is the fine structure 
constant), and 2HB is the damping term related to the 
cosmological expansion. The validity of our approxima- 
tion can be easily verified by noting that the magnetic 
diffusion length Ldiff(t) = •\/i/47rK exceeds safely the 
characteristic length of the relevant process, namely the 
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neutrino mean- free-path i^{t) '-^ 1 s {t/td)^^"^. Using 
the expression of the electric conductivity for relativis- 
tic plasma presented above we find that the condition 
Ldis{t) ^ ^iy{t) is comfortably fulfilled at any time es- 
sential for the presented mechanism. We also checked 
that relativistic corrections to Eq. (Q) do not affect signif- 
icantly our results even if u ~ 1. Thanks to the absence of 
magnetic diffusion, the first term in the r.h.s. of Eq.(|^) 
ensures an exponential amplification of any preexistent 
seed field which is accounted by the the term V x Jext 
in Eq.( Interestingly, neutrino diffusion in the pres- 
ence of a inhomogeneous neutrino-antineutrino asymme- 
try naturally provides such a seed. The relevant point 
here is that, due to charge symmetry breaking in the 
standard model, the neutrino-electron and the neutrino- 
positron weak cross sections are different. As a conse- 
quence electrons and positrons are subject to different 
forces in the presence of a net flux of neutrinos and this 
gives rise to a non-zero electric current. The latter can 
be estimated as follows. The equation of motion of elec- 
tron (or positron) in plasma is Pe — '^uO've^P — Pe/Te 
where = 5nv is the neutrino flux with Sui, being 
an excess of neutrino number density, Ap ^ T is the 
transferred momentum, and Te = {aTn^)~^ is the char- 
acteristic damping time due to interactions with pho- 
tons, (Tt = 87ra^/3mg is the Thomson cross-section and 
— 0.24T'^ is the photon number density. 
Since the ve'' and ve^ cross-sections are different, 
a^^- = G|s/7r « lOGl^T^ and CT^e+ = G\s/in « 
3G^T^ (barring finite electron mass corrections), the 
drift velocities of electrons and positrons would be dif- 
ferent and a non-zero electric current would be induced 
by neutrino flux. Using the previous expressions the 
difference of (e~ — e+)-velocities can be estimated as 
Au « 2 • 10-'^^{5n^/n^){T/MeVf . Correspondingly the 
electric current induced by neutrinos is given by: 

which corresponds to a seed field with strength B^^^*^ w 

lO^^'^ (iv^)^ ^ t/A ^ 1. We are now able to solve 
numerically Eq.(^) by coupling this equation to Eqs.(|[^. 
Instead of the absolute value of B it is convenient to fo- 
cus on the ratio b = \eB\/T'^ which is a constant quantity 
for a frozen-in MP in the absence of entropy production. 
From Fig. 1 it is evident at a glance the huge amplification 
of b which takes place in the short interval of time during 
which fluid "turbulence" is active. A suitable cutoff has 
been used in our computation to account for the satu- 
ration of the amplification process which has to come in 
when the MF approaches energy equipartition with the 
fluid motion, i.e. when B^/Ai: ^ {p + p)v^. We think 
that a more elaborate treatment of non-linear MHD ef- 
fects might only slightly change the knee of b{x) without 
affecting significantly our main results. The final value 



of b depends on the rapidity of the amplification process. 
Since b{x) cx exp {v{x) x/\), it is clear that equipartition 
will be reached more quickly, i.e. when v is yet not sup- 
pressed by viscosity, for fluctuations having large ampli- 
tude and small sizes. In Fig. 2 we present the final value 
of 6 as a function of for several values of [Sn^/riu). 
In principle the requirement of a successful BBN put an 
upper limit on the strength of the MFs, which is roughly 
b ^ 0.1 |9|. This bound, however, was obtained under 
the assumption that the MF was uniform over the Hub- 
ble volume at BBN time and, depending on the details 
of the MF power spectrum, it may not apply to our case. 
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FIG. 2. The final ratio b = \eB\/T^ on the co- 
moving length scale Ad for three values of the neutrino 
number density contrast: starting from the right curve 
i5n,/n^)x^ = 1, 10-\ 10"^ 

In principle inhomogeneous MFs can be subject to dis- 
sipation due to their back-reaction on the plasma. In 
fact, tangled MFs give rise to an anisotropic pressure that 
could induce plasma oscillations (mainly Alfven waves) 
about a force-free configuration. These oscillations might 
be rapidly damped because of the high fluid viscosity 
which implies a dissipation of magnetic energy ||l0|] . In 
our case, it is easy to verify that during the neutrino de- 
coupling process the oscillation frequency of Alfven waves 
is always smaller than the characteristic rate at which 
MFs are produced which is < H so that these MHD 
modes are not excited. MF dissipation could also take 
at photon decoupling. However, it was showed in Ref. 
|]lO| that such an effect is negligible for MFs extending 
on galactic scale if their present time strength is smaller 
than 10^^ Gauss. 

The generation mechanism that we have discussed 
gives rise to MFs of intensity Bq = b 8 x 10~^ Gauss 
at the present time with a coherence length Aq = 
Arf rnitd) (Td/To) ~ 10^ A^ pc. Galactic MFs are ob- 
served with characteristic strength of the order of 1 
extending over scales ~ 1 kpc ||l|. Taking into account 
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flux conservation during the protogalaxy collapse, the 
primordial origin of galactic fields would require a proto- 
galactic field with the strength ~ 10~^° Gauss and the 
coherence length of 0.1 Mpc [0J|]. Although this scale 
is much larger than the coherence length predicted by 
our model, it is natural to expect that some homogeniza- 
tion could take place during galaxy formation. Since the 
field orientation is random over scales larger than Aq, the 
predicted mean field on the protogalactic scale will be 
obtained by a suitable volume average M 



5(0.1 Mpc) ~ Bq 



Ao 



0.1 Mpc 



3/2 



10~i°6 A//2 G. (6) 



From this equation and our previous results (see Fig. 2) 
we find that galactic MFs may be a product of neutrino 
number fluctuations with the amplitude ^ 1 extending 
over scales comparable to the Hubble horizon at neutrino 
decoupling. It is remarkable that MFs with this inten- 
sity may have produced observable effects on the CMB 
and have interesting consequences for structure forma- 
tion . Fluctuations with a smaller amplitude could still 
have played a role in the generation of galactic MFs if a 
galactic dynamo, or even much less efficient amplification 
processes, took place during or after galaxy formation. 

Another interesting consequence of stirring the pri- 
mordial plasma by neutrino inhomogeneous diffusion is 
the production of gravitational waves (GW). The gen- 
eration of a cosmic background of GW by turbulence 
produced at the end of a first-order phase transition 
was discussed in Ref. [|ll|. In our case GW with the 
largest amplitude are produced with a comoving wave- 
length Ao = Xd 10^ pc which corresponds to the present 
time frequency t^o — 10^^^^"'^ Hz. The GW production 
took place when i'iy(x*) = A(x*), i.e. = ^/Xd, and it 
lasted for a time interval comparable to iJ(a;*). Follow- 
ing ref. 1 1 1 1 we estimate the energy density parameter of 
GW with frequency ujo to be 



10-^ H{x,) X{x^)v^{x^) 
3 X IQ-'^Xl^^ v'^ix,) . 



(7) 



From our previous results it follows that a GW back- 
ground produced by neutrino number fiuctuations of am- 
plitude < 1 with Xd < 10^^ may be detectable by future 
GW space based observatories . As tangled MFs can 
also act as a source of GW |l^, a further, and perhaps 
dominant, contribution to the GW background may come 
from MFs produced during the decoupling process. 

Wc conclude this Letter by observing that large neu- 
trino number density fluctuations, as those required to 
power the effects discussed above, might be generated 
by several mechanisms. One possible way to gener- 
ate inhomogeneous and large lepton asymmetry at large 
scales together with a small baryon asymmetry could be 
achieved ||lj,|l^ in the frameworks of Affleck-Dine fl^ ] 



baryo(lepto)-genesis scenario. Another very interest- 
ing mechanism based on active-sterile neutrino oscilla- 
tions was proposed recently by DiBari |]l^. Such a 
mechanism naturally gives rise to domains where neu- 
trinos, or antineutrinos, of a given species are strongly 
converted into sterile neutrinos (or sterile antineutri- 
nos), hence dnv/rii, ~ 1. The typical domain size 
is determined by the neutrino mean free path at the 
critical temperature at which neutrino conversion takes 
place ~ 15 MeV {\5m'^\/eW^f^ , which implies Ad ~ 

10^'' (|(Sm^|/eV^) A nice feature of this scenario is 
that it practically does not need to invoke new physics, 
except for only one that there exists a sterile neutrino 
mixed with an active one. In this case very small inhomo- 
geneities in the baryon number density, which is known 
to exist in the early universe, would give rise to a very 
strong amplification of initially negligible lepton asym- 
metry. Since the initial excess of energy density of active 
neutrinos is exactly compensated by the deficit in sterile 
ones such a model does not suffer from a possible dis- 
tortion of approximate isotropy of the cosmic microwave 
background radiation. Both mechanisms discussed here 
produce isocurvature fluctuations in the neutrino fluid. 
Chaotic vector and tensor perturbations, which may be 
absent initially, are produced when neutrino start to dif- 
fuse according to the mechanism discussed in this Letter. 
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